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We present improved circuits for the control-control-phase (Toffoli) gate and the control-swap (Fredkin) gate 
using three and four global two-qubit gates, respectively. This is a nearly double speed-up compared to the 
conventional circuits, which require five (for Toffoli) and seven (for Fredkin) conditional two-qubit gates. We 
apply the same approach to construct the conditional four-qubit phase gate by seven global two-qubit gates. 
We also present construction of the Toffoli and the Fredkin gates with five nearest-neighbour interactions. Our 
constructions do not employ ancilla qubits or ancilla internal states and are particularly well suited for ion qubits 
and for circuit QED systems, where the entangling operations can be implemented by global addressing. 


PACS numbers: 03.67.Lx, 03.67.Ac, 37.10.Ty, 32.80.Qk 

I. INTRODUCTION 

The three-qubit Toffoli and Fredkin gates are arguably the 
most important highly-conditional quantum gates in quantum 
information science because they are the key enabling ingre¬ 
dient in quantum error correction mm. It is known that 
these gates are universal: any quantum computation can be 
constructed entirely by Fredkin or Toffoli gates, combined 
with the one-qubit Hadamard and phase gates 13 II. The 
Toffoli gate or the closely related control-control-phase gate 
(cc-phase gate) has been experimentally demonstrated with 
nuclear magnetic resonance (NMR) 0, linear optics 00, 
trapped ions 0, and recently with superconducting qubits 18- 
HOl . The Fredkin gate has been demonstrated with NMR E). 

The standard decomposition of the Toffoli gate in the circuit 
model of quantum computation uses six CNOT (or c-phase) 
gates sum It can be constructed also with five conditional 
two-qubit gates — two CNOT gates and three controlled-V 
gates, where V is the \/NOT gate | l_2. !_3| (see Appendix | a|. 
The Fredkin gate can be obtained from the Toffoli gate if sur¬ 
rounded with two CNOT gates Recently, it was shown 
that five conditional two-qubit gates are not only sufficient but 
also necessary for the implementation of the Toffoli gate m 
and the Fredkin gate Ifl4l[l6 |. 

The speed of quantum circuits is usually determined by the 
number of conditional gates and the objective is to minimize 
this number. Over the last years, several alternative methods 
for more efficient construction of the three-qubit gates have 
been proposed and demonstrated. Following a theoretical pro¬ 
posal by Ralph et al. na. Lanyon et al. 0 constructed the 
Toffoli gate in a linear optics experiment with only three two- 
qubit gates by using a qutrit target, i.e., with an ancilla state 
added to the target qubit. Ivanov and Vitanov lfl8l proposed 
a method for constructing highly-conditional C"-NOT gates, 
which uses composite pulses, without two-qubit gates. Monz 
et al. 0 have demonstrated experimentally the Toffoli gate 
with calcium ions in a linear Paul trap by using a sequence 
of 15 laser pulses; 13 of these pulses were on the first blue- 
sideband transition and they performed reversible mapping 
between electronic and vibrational motion, as in the Cirac- 
Zoller’s CNOT gate OWTIl . 

In this paper, we propose realizations of the three- and the 
four-qubit conditional-phase gates and the Fredkin gate using 


global entangling gates and no ancilla internal states or ancilla 
qubits. Two types of physical systems are considered for the 
realization of the entangling gates, where i) all qubits interact 
simultaneously (global interaction), or where ii) each qubit 
interacts with its nearest neighbours only (nearest-neighbour 
interaction). The three- and the four-qubit phase gates are con¬ 
structed with only three and seven global interactions, respec¬ 
tively, and the Fredkin gate - with four. Each of the global 
two-qubit gates is implemented in a single step, by global ad¬ 
dressing of all qubits, which amounts to simultaneous two- 
qubit phase gates upon each pair of qubits. These realizations 
are particularly well suited for trapped ions in a laser-driven 
linear Paul trap 1(2214241 . and in a microwave-driven trap with 
magnetic-field gradient 12514281 . We note that this speed¬ 
up does not contradict the earlier findings that five two-qubit 
gates are both necessary and sufficient for the three-qubit Tof¬ 
foli and Fredkin gates □SHE!, because each global gate used 
here embraces three two-qubit gates. The total number of two- 
qubit gates is therefore larger than five; however, they are con¬ 
veniently grouped into fewer global gates. 

With nearest-neighbour interactions the three-qubit c-phase 
and Fredkin gates require five entangling gates. These realiza¬ 
tions are suitable, e.g., for circuit cavity QED 1(291 1. where the 
nearest-neighbour coupling arises naturally from the hopping 
of photons between adjacent cavities in the regime of photon 
blockade l30l . 

Global addressing is the key tool in the Sprensen-Mplmer 
two-qubit c-phase gate ED and in their proposal for genera¬ 
tion of GHZ states |[32l . It has been used recently in the con¬ 
struction of the Toffoli gate by only three global conditional 
operations l33l . Here we extend this method to other gates 
and also present constructions with global addressing in the 
important case of nearest-neighbor interaction. 

Physical realizations with trapped ions are also discussed. 
In the laser-driven Paul trap we use bichromatic laser fields, 
the two components of which are tuned at a certain detuning 
from the first red and blue sidebands of the qubit transition. 
These fields address simultaneously all qubits at each step. 
With suitably chosen laser phases these bichromatic fields 
create effective spin-spin interactions between the qubits and 
cancel the dependence on the phonon number. The thereby 
designed Toffoli and Fredkin gates are thus both faster and ap¬ 
plicable to ions in a thermal state of motion, like the Sprensen- 
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M0lmer c-phase gate ED. 

In the microwave-driven magnetic-gradient trap I [2514281 . 
the gates are implemented in a particularly simple fashion be¬ 
cause the required spin-spin couplings arise naturally from the 
magnetic-field gradient. The conditional two-qubit gates are 
implemented merely by letting the system evolve freely for a 
certain time duration. 


II. FREDKIN, TOFFOLI AND THREE-QUBIT 
C-PHASE GATES: MATHEMATICAL CONSTRUCTION 


state 11), but are left unchanged otherwise. In the basis |T]i the 
Fredkin gate has the matrix form 

' 1 0 0 0 0 0 0 0 ' 

0 1 0 0 0 0 0 0 

00100000 
00010000 
T = 00001000 ■ 
00000010 
0 0 0 0 0 1 0 0 

0 0 0 0 0 0 0 1 


In the Toffoli gate, there are two control qubits and a target 
qubit: the target qubit is inverted if the two control qubits are 
in state |1), but it is left unchanged otherwise. In the three- 
qubit basis. 


41000), 1001), 1010), 1011), |100),|101),|110),|111)}, (1) 


the Toffoli gate has the matrix form 
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( 2 ) 


Here we introduce improved quantum circuits, which pro¬ 
duce the three-qubit cc-phase gate, the four-qubit ccc-phase 
gate and the Fredkin gate with a minimal number of condi¬ 
tional operations. We consider circuits based on two types of 
interactions: global interactions, where all pairs of qubits in¬ 
teract collectively, and nearest-neighbour interactions, where 
each qubit interacts with its nearest neighbours only. Both re¬ 
alizations can be produced by fields addressing uniformly all 
qubits. 


A. Implementation with global interactions 

We begin with the design of circuits where the conditional 
gates are produced by global interactions, which we call G- 
gates. 


and the three-qubit conditional phase (cc-phase) gate is 


1. Toffoli and cc-phase gates 


1 0 0 0 0 0 0 0 

0 1 0 0 0 0 0 0 

0 0 1 0 0 0 0 0 

0 0 0 1 0 0 0 0 

0 0 0 0 1 0 0 0 

0 0 0 0 0 1 0 0 

0 0 0 0 0 0 1 0 

0000000 -1 


(3) 


By using the method outlined in Appendix [B] we have 
constructed a quantum circuit of the three-qubit conditional 
phase gate with three G-gates only, which can be verified by 
direct inspection. The circuit is shown in Fig. IIA 1 (top). 


Besides the G-gates, it consists of n/2 pulses (boxes), and 
phase gates (ovals). See Appendix [A] for the complete list of 
notations. In operator form (to be read from left to right), the 
circuit reads 


Because the cc-phase gate is symmetric, and all three qubits 
have the same role (no distinction of target and control gates) 
hereafter we shall be concerned only with it. As it is well 
known, the Toffoli gate can be obtained from the cc-phase 
gate with two Hadamard gates. 



HZ3-H1F 


Given the cc-phase gate, any of the three qubits can be 
turned into a target qubit by the application of the two 
Hadamard gates on it. Therefore, the implementations of the 
cc-phase gate below can be directly translated into Toffoli 
gates. 

In the Fredkin gate, there is one control qubit and two target 
qubits: the target qubits are swapped if the control qubit is in 


BiG(f)B 1 <l> l (f)B I G(f)B 1 G(f)<l> 1 (f)0 2 (f)03(f), (4) 

where 


to* 

=Jh* 

T 

II 

(5a) 


(5b) 

G(0) = e i *(°f 0 2+ 0 2 0 3 + 0 l 0 3). 

(5c) 


This realization demands only three conditional gates, four 
phase gates (in two phase groups, hence the T-depth ll34l is 
2) and four n/2 rotations. Because the gate implementation 
time is determined primarily by the number of conditional 
two-qubit gates this implementation is faster by a factor of 2 
than the traditional implementation with six CNOT gates (see 
Appendix [A). 

We note that a similar circuit for the Toffoli gate with three 
global conditional operations has been proposed by Ref. ]33l 
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where Qk represents a it/A rotation on qubit k, Q k = e ' * °i-. 


3. Four-qubit ccc-phase gate 

The same approach can be used to constmct efficient cir¬ 
cuits for higher-conditional gates. Here we present an exam¬ 
ple for the four-qubit ccc-phase gate, which in operator form 
reads 



= 1 — 2 | 1111 )( 1111 | . 


Figure 1: Global realization of the three-qubit c-phase gate (top) and 
the Fredkin gate (bottom). 


who have used five global single-qubit n/2 and it/A rotations 
and three local phase gates rather than the four local n/2 ro¬ 
tations and the four local phase gates used here. 

Due to the commutation of the operators of the G-gate can 
be represented as a product of three O z a 7 k factors, G(<j>) = 
e i| l ICT iO'2e^ <I 2 tT ! e ^ CT r (I 3. Because of the property 


j2 _ , 1 Y‘ 


N aV 

j<k °J °k > 


( 6 ) 


where Jp = a j (P = x, y,z) is the collective spin pro¬ 

jection operator, the G-gate is equivalent, up to a phase factor, 
to a similar gate proposed by Mplmer and Sprensen (32l for 
creation of GHZ states (35). Such global gates can be im¬ 
plemented by a single driving field interacting uniformly with 
all qubits of the desired target gate. The Jr gate is a part of 
a universal set of gates proposed by Nebendahl et al. (33) as 
an efficient alternative of the standard toolbox involving the 
CNOT or the two-qubit c-phase gate. 

Rather than using rotations and phase gates, all circuits 
from this paper can be written with Hadamard gates and T- 
gates, traditionally used in information theory. This can be 
achieved by using the relations (A2| >. For the cc-phase gate 
we have (in operator form, to be read from left to right): 


Hi G( f )Hi 7j*tf! G( f )Hi G( §) 7j T 2 T 3 , (7) 


This gate can be obtained with the following sequence of in¬ 
teractions, containing 7 global G-gates (to be read from left to 
right): 


PG(f)4>(f)H 123 G(f)<I>(-f)P4> 4 (-^)G(f)PG(f)x 

x® 4 (-f)W(|)®(^)TO(f)®H)»4(f)x (9) 

xG(f)*(-?)/*G($)* 4 (f)*(?f), 

where 


p 12 3 = <r‘?K+°2 +0 3), 

(10a) 

p _ e -if (of+t^+of+oj) ) 

(10b) 

<J>(0) = e “' < H CT f+ 0 2 + 0 3+ cr 4) ; 

(10c) 

G(0) = g ! W CT i cr 2+ CT f CT 3+ CT i CT 4+®2 CT 3+ cr 2 <T 4+ <T 3 0 4). 

(10d) 


Numerical simulations suggest that 7 is the minimum number 
of global gates G(0). If we allow for a single non-global two- 
qubit gate, this number is reduced to 6. Barenco et al. pro¬ 
posed a scheme with 13 conditional gates - 6 CNOTs and 
7 c-phase gates 02, which is twice as many as we propose. 
Other proposals exist, which prescribe fewer conditional gates 
m . However, they use either ancilla qubits or ancilla states. 


B. Implementation using nearest-neighbour interaction 


where the index indicates the qubit number. 


If the qubits interact only with their nearest neighbours, as 
in circuit QED, a global addressing produces the gate 


2. The Fredkin gate 


N((j>) = e'0( CT l O 2+ O 2 O 3) ) 


(ID 


We have constructed a circuit for the Fredkin gate <(4)) with 
four global two-qubit gates, shown in Fig. IIA 1 1 (bottom). 
Here qubits 2 and 3 are swapped conditionally on the state of 
qubit 1. In operator form (to be read from left to right): 


F 3 03 (f)G(f)P 2 F 3 ‘l>i(f)‘I> 2 (f)‘I> 3 (f)G(f) J P 2 G(f)x 
x0 2 (f)2 2J P 3 <I> 3 (f)G(f) J P 3 <I> 3 (f), (8) 


rather than the G-gate used above. Here the middle qubit 2 
interacts with the outer qubits 1 and 3, while qubits 1 and 3 do 
not interact directly. The /V-gate produces less entanglement 
than the G-gate; hence the quantum circuits with A'-gates are 
longer than with G-gates. We point out that the A'-gate, like 
the G-gate, can be achieved in a single interaction step by driv¬ 
ing all qubits globally with the same field. 
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Figure 2: Nearest-neighbour realization of the three-qubit cc-phase gate (top) and the Fredkin gate (bottom). 


1. The cc-phase gate 


The cc-phase gate can be obtained using the circuit, shown 


in Fig. IIB 1 (top). Now five /V-gates are necessary, as op¬ 


posed to the three G-gates used by the global implementation. 
We note here that the number of conditional gates is the same 
as in the standard circuit model 141 1121 : however, the standard 
circuit (see Appendix [A| uses conditional two-qubit gates be¬ 
tween each pair of qubits, rather than between nearest neigh¬ 
bors only. Here we use only nearest-neighbors interactions. 


2. The Fredkin gate 


We have also constructed a circuit for the Fredkin gate. 


given in Fig. IIB 1 (bottom). The control qubit is 1 and the 
target qubits to be swapped are 2 and 3. Like the cc-phase 
gate, the Fredkin gate demands five conditional gates. 


III. PHYSICAL REALIZATION WITH TRAPPED IONS 


Below we consider a few natural physical realizations of 
the G-operator with trapped ions. Two architectures are con¬ 
sidered - laser-driven traps and magnetic-gradient traps. Be¬ 
cause the G-operator bears the highest cost in time duration 
we focus our attention at its construction. It is noteworthy that 
the G-operator is symmetric upon permutation of the indices 
of the <7 Z operators, and it involves all qubits. The important 
implication from here is that, in laser-driven traps, it can be 
implemented in a single shot by a bichromatic pulse, i.e. by 
addressing all three ions simultaneously with the same field, 
in a global fashion. In magnetic-gradient traps, it is delivered 
simply by free evolution. 


A. Realizations with bichromatic laser interaction 


Z axis. We assume that the ion crystal interacts with a bichro¬ 
matic laser field along the transverse x direction with laser fre¬ 
quencies C0r = cqq — (Oc.m. + 8 (“red”) and cob = ®o + © c . m . — 8 
(“blue”) tuned near the center-of-mass (c.m.) mode with 
detuning 8, (ft) cm |5|). The interaction Hamiltonian in the 
Lamb-Dicke limit (rj <C 1) and the rotating-wave approxima¬ 
tion is given by f37l (h= 1) 

Hi = EL +ae - iSt+i *). (12) 

Here gk{t) = /\/N is the individual time-dependent 

spin-phonon coupling for ion k, tj is the single-ion Lamb- 
Dicke parameter, and £2*(f) is the Rabi frequency “unit”, a' 
and a are creation and annihilation operators of the c.m. vibra¬ 
tional mode and = e~ ul> k o k with cr^ being 

the spin raising and lowering operators for ion k. In terms of 
the blue and the red laser phases <pj? and <p k , the spin and mo¬ 
tional phases in Eq. ( fl2| ) are defined as <p^~ = \((p k + cpji) and 
cp { 7 = \(<p k — (p[,), respectively. Hereafter we assume that the 
blue and the red lasers are counterpropagating, which implies 
that (p k = 0 and (p k = <jsuch that ) = a k . 


1. The a x gate 


Let us consider the case of N ions addressed simultaneously 
with a suitably detuned bichromatic field, as in the S 0re risen- 
Mqlmer two-qubit c-phase gate ED. Setting g = g k (t) and 
<p = <j \ in Eq. m the propagator becomes 


U = D(a)exp 



(8t — sin 8t)J^ , 


(13) 


where D(a ) = e aa “ +a is the displacement operator with a = 
2gJ x e~ l ^(\ — e ,St )/8 being the displacement parameter. At 
time 


We consider a linear ion crystal of N identical two-state ions 
with transition frequency coo confined in a Paul trap along the 


_ 2 K 
Igate — “g - 7 


( 14 ) 
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the motional displacement in Eq. ( fl3| vanishes, a = 0, while 
the spin component acquires a nonzero geometric phase shift, 

U =exp^' 8 ||-i^ . (15) 


Using the property ([ 6 ]) we finally obtain 


U = exp 



(16) 


Note that in all formulas above, N = 3 for the three-qubit gates 
and N = 4 for the four-qubit gates. Apart from the unimpor¬ 
tant phase factor, the propagator of Eq. ( p~6] > is exactly the 
G-gate of Eq. ( |5cj ) needed in our implementations, with the 
relation 


0 


_ 4ng 2 

= 


(17) 


We note here that this implementation of the G-gate with 
the o x couplings is equivalent to an implementation in the o z 
basis (described in the charts above) rotated at an angle of n/4 
(see details in the Appendix). 


2. The d z gate 

The alternative scheme is based on collective Raman type 
interaction with frequency difference close to the vibrational 
frequency rather that the qubit frequency. We assume that the 
linear crystal of N ions is uniformly addressed with two Ra¬ 
man beams with wave vector difference AA pointing along the 
transverse x direction and laser frequency difference Aft)| = 
®c.m. — S. The Hamiltonian describing the laser-ion interac¬ 
tion in the Lamb-Dicke limit is given by 1381139 1 

H { = gJ z (a f e iSl +ae~ iSt ). (18) 

The unitary propagator corresponding to the Hamiltonian ( p~8| ) 
is identical in form to (f]~5j) by replacing J x —> J-. 


B. Realization with magnetic-gradient ion traps 

In this ion-trap scheme a magnetic field is applied such that 
the ions experience a field gradient along the chain. This cre¬ 
ates a magnetic gradient induced coupling (MAGIC), which 
makes microwave or rf radiation effective for addressing and 
coupling internal qubit states to the vibrational motion of the 
ions. Thus the magnetic-gradient ion trap does not require 
laser light for single- and multiple-qubit quantum gates l25l . 
The Hadamard gate is achieved with a microwave or rf k/ 2- 
pulse resonant with the qubit transition frequency, while the 
phase gate can be obtained by choosing the phase of the driv¬ 
ing field. The conditional gate ( |5c| can be obtained in a very 
simple fashion, as described in the rest of this section. We 
follow the method proposed in Ref. |40| . 


In a frame rotating with the ions’ internal states and to sec¬ 
ond order in the ions vibrational motion the Hamiltonian can 
be written as li25l 

H = Y^,= i - \ ( 19 > 

where a\ and a n are respectively the creation and annihilation 
operators for the nth vibrational mode. The first term is a sum 
of the energies of all N vibrational states. The second term 
represents long-range pairwise spin-spin coupling, whereby 
spin j is coupled to spin k with the coupling coefficient ll26l 

Jjk= (ge ^ )2 (A~V (20) 

Here A is the Hessian of the trap potential taken at the equilib¬ 
rium positions of the ions, b is the amplitude of the magnetic 
field gradient, g e and Lip, are respectively the electron g-factor 
and Bohr’s magneton. This spin-spin coupling is only weakly 
sensitive to thermal excitation of the ion string. Thus, usually 
just Doppler cooling is sufficient to avoid unwanted thermal 
effects on the coupling constants Jp, HTl . 

The second term in Eq. m which describes the spin 
states, is decoupled from the first term, so the spin states can 
be described independently from the motional states. The 
propagator corresponding to this spin term is 

U( T) = exp (^Ylj<k=i J jk a j a i) ■ ( 21 ) 

It describes the free evolution of the N spins (with no mi¬ 
crowave or rf driving), which are only exposed to the magnetic 
field gradient for time period f. 


1. Equal couplings 

For equal couplings, = /, the propagator of Eq. CCQ> 
reduces to the conditional gate G(0) at time T = 20 /J, 

G(0) =U{2<j>/J). (22) 

Thus the implementation of G(0) requires no external driving. 
It simply takes place as a result of the free evolution of the 
spins for time 20 /J. Hence, in the magnetic-gradient traps the 
Toffoli gate and the Fredkin gate can be realized as sequences 
of microwave or rf pulses (producing the single-qubit gates), 
applied at certain times BOl . 


2. Unequal couplings 


Our method can be applied for different couplings Jp, as 
well. In the linear Paul trap with harmonic potential, for 
example, where J \ 2 = J\ t / Ji 3 , the three-qubit conditional 


phase gate is obtained with the circuit, shown in Fig. IIIB2 
In operator form (to be read from left to right): 


P 1 M(f)P 1 T> l (f)P 1 (0 1 )M(0 1 )P 1 (0 2 )M(0 2 )x (23) 

xPi(0 3 )Oi(f)<I> 2 (f )4> 3 (f). 
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Figure 3: Realization of the three-qubit cc-phase gate with unequal 
couplings J For uniform J j k this circuit reproduces the circuit from 
Fig. fTTXU ftop). 


Now the conditional gate ( pc) is defined as M ( (j)) = 
U(2<j>/J l2 ) and P k {6) = e~^ a i. We have fa = 0.375?:, fa = 
0.258?:, 0i = 1.0327:, 0 2 = 0.484?:, 0 3 = 1.484?:. 

Note that we have relabelled the ions from 1,2,3 to 2,1,3. 
Thereby, for uniform Jj k , the above circuit is reduced to the 
circuit from Fig. IIA l| (top), as expected. 

The corresponding circuit for the Fredkin gate is not given 
as it is more cumbersome. 


IV. DISCUSSION 

We note here that several previous works have paid special 
attention to the single-qubit phase gates involved in the con¬ 
struction of the Toffoli gate because in some quantum com¬ 
puting platforms these gates are the ones that are most diffi¬ 
cult to construct. As mentioned above, the standard imple¬ 
mentation of the Toffoli gate uses six such gates in five groups 
(meaning that two of them can be applied simultaneously to 
different qubits) mm. These numbers have been reduced 
to four phase gates in two groups by using ancilla states and 
six CNOT gates. Amy et al. l34l reduced the T-depth to 3. In 
all of our implementations proposed here we have the same 
numbers — four single-qubit phase gates in two groups — 
but without ancilla states and using 3 global two-qubit gates. 
Hence in the information theory language, our construction of 
the Toffoli gate has T-depth equal to 2. Therefore our con¬ 
struction of the Toffoli gate is as compact as the best T-depth 
of achieved previously. However, we emphasize that the T- 
depth is of little interest to implementations with trapped ions 
because for them the number of two-qubit conditional gates 
(Clifford gates in information science language) is of primary 
concern because this number determines the speed of the total 
gate. 


V. CONCLUSIONS 


2 faster than conventional proposals with six CNOT gates (see 
Appendix [A). The Fredkin gate uses four global conditional 
gates. The four-qubit ccc-phase gate uses seven global condi¬ 
tional gates or six, if one of the conditional gates is not global. 
Proposals exist, which prescribe fewer such gates, however 
resorting to ancilla states or ancilla qubits. We point out that 
no ancillas are used in our circuits. 

We also give various physical realizations. Global condi¬ 
tional gates are perfectly suited for systems such as ion traps, 
where ions interact equally with an applied external field. Two 
types of traps are considered. In the standard laser-driven 
traps, each conditional gate is implemented using a single 
bichromatic pulse; two particular realizations are given. In 
the microwave-driven traps the conditional gates are achieved 
simply from the/ree evolution, i.e. with no external field driv¬ 
ing - we only let the ions evolve freely for a predetermined 
time period under the influence of the magnetic-gradient field. 

Nearest-neighbour interaction is suitable, e.g., for cirquit 
QED, where the nearest neighbour coupling arises naturally 
from the hopping of photons between adjacent cavities in the 
regime of photon blockade {30). Our circuits of the Toffoli 
and the Fredkin gate prescribe five nearest-neighbour condi¬ 
tional gates. We note that all such gates are also global in the 
sense that they can be achieved by driving all qubits simulta¬ 
neously; therefore they can be realized in a single interaction 
step, too. 

Finally, we point out that our numerical methods cannot 
give a proof that our circuits are optimal. But we note that 
after extensive search no better circuits using the described 
gate set were found. 
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Appendix A: Notation for gates 


The traditional implementation of the Toffoli gate uses 6 
CNOT gates and 10 single-qubit gates PflfTZI, 


- 0 - 




404^ 


In this work we have proposed mathematical constructions 
of the Toffoli gate (or the closely related two-qubit cc-phase 
gate), the four-qubit ccc-phase gate and the Fredkin gate us¬ 
ing global two-qubit gates. The global gates employ either 
i) global interaction, where all pairs of qubits interact collec¬ 
tively, or ii) nearest-neighbour interaction, where each qubit 
interacts with its nearest neighbour qubits only. The Toffoli 
gate demands three global two-qubit gates, which is a factor of 


The three-qubit c-phase gate is constructed with the same cir¬ 
cuit but without the Hadamard gates. 





r t 4r i '4j^ 
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The traditional single-qubit gates are 


-Sr- 


'or 

1 0 

Pauli O x 

(Ala) 



0 -i~ 
i 0 

Pauli a y 

(Alb) 



1 0 

0 -1 

Pauli cr 

(Ale) 

-10 

1 

V2 

' 1 1 

1 -1 

Hadamard 

(Aid) 



'10' 

0 i 

n/A phase gate 

(Ale) 



1 0 

0 e in ! A 

7t/ 8 phase gate 

(A If) 


IIe iaa H = e~ ia(rX 


cos a — tsina 
—tsina cos a 


(A7) 


(Hi ®H 2 )e~ iaa i®°2 (Hi ®H 2 ) = e _ ‘ aCT i ®°2 (A8) 

(Hi <g> H 2 )e~ iaa (H i ®H 2 ) = e ~ iaa ^ (A9) 

Appendix B: Numerical procedure 

We follow the procedure as described in Ref. |40|| . First, 
we form a generic circuit containing Nq two-qubit gates B: 

F = B(4> Ng )L N g ... B(fo)L 3 B(h)L2B(fa)Li, (Bl) 


-0- 


J_ 1 1 
Vi i 1 


\/NOT gate 


(Alg) 


In a closed qubit, all unitary transformations are SU(2), i.e. 
with unit determinant. In this case, we use the physically rel¬ 
evant gates 




050 


0 — i 
—i 0 


l 

V2 


1 — i 
—i 1 


0 

0 


K pulse, e ? 

n /2 pulse, P = e~ l 2 3 4 ^ a * 

phase gate, <f>(<f>) = 


where implements a generic local operation to each qubit, 
Lk G S U (2) (g> SU (2) Cg) S U (2), and B represents one of the con¬ 
ditional gates G, N or M. An element h from SU(2) can be 
constructed, up to an unimportant global phase, in the follow¬ 
ing way: 


h = <P((pi)ex V (-tOG x )<t>((p 2 ). (B2) 

This represents a pulse with area 0, surrounded by two phase 
gates. Because the phase gates <t>(<p) commute with the oper¬ 
ator B, Of (pi) from L* can be combined with <P((fh) from Lj ,, | 
into a single phase gate. Therefore, without loss of generality, 
we take all (p 2 to be 0. 

Then we proceed with a numerical minimization of the dis¬ 
tance 


Unlike the preceding gates, these gates can be implemented 
in a closed qubit, without ancilla states. Obviously, 


O x = ie-‘2°\ 

(A2) 

S = e in ^(n/ 4), 

(A3) 

T = e iK/s ®(ji/8), 

(A4) 

H = ;'T>( 7r/4)e“ ! ? CT '4>(a:/4). 

(A5) 


Useful relations between x and z bases: 


He- iaa 'H = e~ ia& 


~ ia 0 
0 e ia 


( A6) 


D = L\ F ij~Fij |> ( B3 > 

hi 


with F being the target gate. We use Newton’s gradient-based 
method to determine the variables ©/., (j)/ and 0/ ( , which yield 
D = 0. From all solutions obtained, we choose those with the 
fewer number of single-qubit gates. To minimize the number 
of the two-qubit gates B((j>), we start from a small number 
/Vg, which we gradually increase, until we reach a solution 
to D = 0. Because we use a local optimization algorithm, 
we iteratively pick the initial values of the variables using a 
Monte-Carlo scheme. 
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